The target problem with evanescent subdiffusive traps 
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We calculate the survival probability of a stationary target in one dimension surrounded by 
diffusive or subdiffusive traps of time-dependent density. The survival probability of a target in 
the presence of traps of constant density is known to go to zero as a stretched exponential whose 
specific power is determined by the exponent that characterizes the motion of the traps. A density 
of traps that grows in time always leads to an asymptotically vanishing survival probability. Trap 
evanescence leads to a survival probability of the target that may be go to zero or to a finite value 
indicating a probability of eternal survival, depending on the way in which the traps disappear with 
time. 

PACS numbers: 82.40.-g, 82.33.-z, 02.50.Ey, 89.75.Da 



O 



o 
o 



> 
in 

o 
o 



O 

o 



X 



I. INTRODUCTION 

The traditional trapping problem involves diffusive 
(Brownian) particles (A) that wander in a medium 
doped with static traps (B) and disappear when th^ 
meet [llHI3,l3. In the traditional iarget problem dUQ, 
on the other hand, one has static A particles and mo- 
bile traps. Both of these problems are described by the 
"reaction" A + B B, but in one case the A's move 
and the B's stand still, while in the other the B's move 
while the ^'s are stationary. Both of these problems 
have a long and active history in the literature. They 
not only represent experimentally observable phenom- 
ena, but they have served as a testbed for theoretical and 
numerical studies and as a starting point for the formula- 
tion of models for more complex systems that have only 
recently been successfully solved analytically. For exam- 
ple, the survival probability of an A particle in a medium 
of B particles when both species are diffusive, first inves- 
tigated numerically in the seminal work of Toussaint and 
Wilczek , was only partially solved analytically 0| 
until the recent full (asymptotic) solution in one dimen- 
sion 0, 0, 01 . These results have also recently 
been generalized to subdiffusive species |l5l| . The survival 
probability of A particles in the reactions A + A ^ A and 
A + A ^ in one dimension when A is mobile is also of 
relatively recent vintage in the history of such analytic 
solutions Hill 113. 

The purpose of this paper is to extend the one- 
dimensional target problem calculations for both diffu- 
sive and subdiffusive traps to the case of traps that them- 
selves disappear in time according to some survival prob- 
ability function of their own (e.g., exponential or power 
law). The decay of the moving traps with time of course 
increases the survival probability of the stationary tar- 
get, and the interesting questions concern the interplay of 
the time dependences of the movement and decay of the 
traps. A related problem was considered in 18], where 
diffusive particles A and traps B and C undergoing the 
explicit reactions (a) A + B ^ B, B + C ^ C, and (b) 



A + B ^ B, B + C were considered using entirely 
different methods. Our methods are equally applicable 
to trap densities that increase with time, but this prob- 
lem is less interesting because it necessarily leads to the 
eventual demise of the target. 

A common characterization of the diffusive motion of 
a particle is through its mean square displacement for 
large t, 
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Here Kj is the (generalized) diffusion constant, and 7 
is the exponent that characterizes normal (7 = 1) or 
anomalous (7 7^ 1) diffusion. In particular, the process 
is sudiffusive when < 7 < 1- Subiffusive processes are 
ubiquitous in nature El IH [H El , and are par- 

ticularly useful for understanding transport in complex 
systems (llisj. 

The problem considered in this paper is a special case 
of a broad class of leaction- subdiffusion processes that 
have been studied over the past decades. One approach 
that has been used to study these processes is based 
on the continuous time random walk (CTRW) theory 
with waiting-time distributions between steps that have 
broad long-time tails and consequently infinite moments, 
ijj{t) ^ t^^^^ for t — > cxo with < 7 < 1. Another 
approach is based on the fractional diffusion equation, 
which describes the evolution of the probability density 
P{x, t) of finding the particle at position x at time t by 
means of the fractional partial differential equation (in 
one dimension) 7. 19-.26..27..2&.29J 
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where is the Riemann-Liouville operator, 
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In this paper we study the one-dimensional target 
problem for a static particle A subject to attack by dif- 
fusive or subdiffusive traps B that may die before reach- 
ing the target A 0. For this pur pose, we general- 
ize the ideas of Bray and Blythe [ij, and of our own 
work ^3 based on a fractional diffusion equation ap- 
proach. While recent work shows that a simple gen- 
eralization of reaction-diffusion to reaction-subdiffusion 
equations in which the reaction and subdiffusion terms 
are assumed to enter additively is not valid in some 
cases [33, this is not a difficulty in our particular ap- 
plication. The difficulties do not arise when the reaction 
process can be translated into a static boundary value 
problem, which is the case for the target (as well as the 
trapping) problem [2^. 



In some cases, asymptotic anomalous diffusion behav- 
ior can be found from corresponding results for normal 
diffusion via the simple replacement of t hy . This can 
be understood from a CTRW perpective because the av- 
erage number of jumps n made by a subdiffusive walker 
up to time t scales as (n) ~ and, in many instances 
the number of jumps is the relevant factor that explains 
the behavior of the system. The simple replacement re- 
sult is evidence of "subordination" (see Sees. 5 and 7.2 
of [^). However, there are other instances where the be- 
havior of subdiffusive systems cannot be found in this 
way. A simple example is the survival probability of sub- 
diffusive particles in the trapping problem (see Sec. 5 
of (6i]). In particular, for systems where competing pro- 
cesses (motion toward target and death) occur according 
to different temporal rules, such a replacement becomes 
ambiguous. This is the case for the problem considered 
here. 

While our analytic results are based on the fractional 
diffusion equation formalism, our numerical simulations 
are based on a CTRW algorithm. These two renditions 
of the problem are expected to differ if trapping events 
are likely in a small number of steps, that is, if the ini- 
tial density of traps is too high. On the other hand, if 
the initial trap density is too low, then the simulations 
to produce valid statistics would take inordinately long 
because trapping events are rare and because the system 
has to be sufficiently large to include many particles. We 
note this as a caveat for our subsequent comparisons. 

In Sec. ^ we present an integral equation for the 
survival probability, which we reduce to quadrature in 
Sec. mil The resulting integral is explicitly evaluated 
for exponentially decaying trap densities (including a 
stretched exponential decay), as well as for trap densi- 
ties that decay as a power law. Not surprisingly, we find 
that a sufficiently rapid decay of the trap density leads 
to a finite asymptotic survival probability of the target. 
Comparisons of our results with numerical simulations 
are also shown in this section. A summary and some 
conclusions are presented in Sec. IIVI 



II. INTEGRAL EQUATION FOR THE 
SURVIVAL PROBABILITY 



We consider a finite interval L containing N = pL mo- 
bile traps B of constant density p initially distributed at 
random, and a single immobile A particle at the origin. 
Following the approach of Bray et al. for diffusive 
traps and our generalization of this approach to the sub- 
diffusive case [l5j, we write the survival probability of 
A as P{t) — exp{— ^o(i)}, where fio{t) is to be deter- 
mined. To find this function, one calculates in two ways 
the probability density to find a B particle at the origin 
at time t, 



dt'po{t')G{t - t') 



(4) 



That the left side is this probability density is obvious. 
On the right side one has the renewal theory expression 
where fio{t')dt' = {—P/P)dt' is the probability that a B 
particle intersected A in the time interval (t' , t' + dt') for 
the first time, and the propagator G{t — t') is the proba- 
bility density for this particular B to be at the origin at 
time t. In one dimension it is given by |3lll32l| 
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(0,1), (1/2,1) 
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where -ff^ '2 is Fox's If-function, whose value at the given 
arguments we have used to write the last equality. In a 
different context than the target problem, Bray et al. 0| 
generalized their approach to a time-dependent density 
p{t) of B. They argue that in place of Eq. Q one now 
has 
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This is the basic equation to be considered in this paper. 

III. THE SURVIVAL PROBABILITY 

To calculate the survival probability of particle A, we 
rewrite Eq. Q explicitly as 
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where we have multiplied both sides by yyiK^, and where 
we have introduced 
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Eq uation © is an equation of Abel of the first kind [s 



fit) = 



r(i-7/2)7o '^^''{t^tr/ 



(10) 



with f{t) = yJ\K^. The solution of this classic equation 
is well known (see Sec. 12 in ,3^ or Eqs. (2.5a) and (2.5b) 
in 0), 



1 
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(11) 



Here, as earlier, qD\ ^^"^ is the Riemann-Liouville frac- 
tional derivative. In our case /(t) — ^ is constant, 
so that 



^ ^ p(t) r(7/2) 



2-1 



It then follows that 



Mo(i) 



r(7/2) Jo 



dzp{z)z 



7/2-1 



(12) 



(13) 



which provides a general solution to our problem for any 
p{t). While Eq. IjlSfl applies to trap densities that grow 
or that decrease or even oscillate in time, the case of 
evanescent traps is the more interesting and the one we 
choose to focus on. 

As a reminder, we note that for traps of a constant 
density p{t) = po the survival probability is given by 



P{t) = exp 



I r(i + 7/2) 



(14) 



As a benchmark, we show in Fig-^a typical comparison 
of this result with simulation results. The agreement is 
clearly good, although a lower initial density run for a 
longer time would lead to even better agreement. 



A. Exponentially decaying trap density 

Suppose that the traps have a finite lifetime r and 
decay exponentially, as in a unimolecular reaction, p(t) — 
Po exp(— t/r). The integral in Eq. IjlSI) immediately leads 
to the solution 



where r(6, x) is an incomplete Gamma function, and 



7U/2 



(16) 




FIG. 1: (Color online) poit) = -In P{t) vs t for non- 
evanescent traps as given in Eq. I|14|l (solid line) and sim- 
ulations (symbols along with error bars). Parameter values 
are 7 = 1/2, po = 0.01, Kj = 1/(2^/^). 



When 7 
duces to 



1, i.e., when the traps are diffusive, this re- 



pait) = ^ipoerf(i 



(17) 



For arbitrary 7 < 1, the survival probability of the tar- 
get in the presence of the subdiffusive traps with finite 
lifetime thus is 



P{t) = exp 



-£~fPo 



1 



r(7/2, t/r) 
r(7/2) 



(18) 



The interesting result here is that the funtion po{t) 
goes to the constant /io(oo) = £-yPo ^ind not to infinity 
as t — > 00. Therefore the survival probability does not 
vanish with increasing time. 



P{t) -> exp 



-iyPO 1 - 



r(7/2)(i/r)i-^/2 



(19) 



We note that is a characteristic distance that mea- 
sures the root mean square displacement of the traps 
during their decay time r. Therefore po{oo) — i-ypo is 
the ratio of this average displacement to the average ini- 
tial distance p^^ between traps. This finite asymptotic 
survival probability, P{oo) = exp(— ^^po)j displays rea- 
sonable qualitative features: it increases with decreasing 
trap lifetime r, and it decreases with increasing initial 
trap density po- That there is a finite asymptotic sur- 
vival probability reflects the fact that if the traps disap- 
pear sufficiently rapidly (which they do if they disappear 
exponentially while the traps move diffusively or subd- 
iffusively), then many traps disappear before they can 
reach the particle, and there is a finite probability that 
the particle remains forever "safe." 

The next two figures show the comparison of simula- 
tion results with our analytic outcome. First, in Fig. [21 
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FIG. 2: (Color online) fj,o{t)/po = — lnP(t)/po vs t for ex- 
ponentially evanescent traps. Solid line: Eq. 1171 . Squares: 
simulation results for a high initial density po — 0.1. X's: 
simulation results for a lower initial density po ~ 0.01. Other 
parameter values are 'y — 1, t — 100, and Ki — D — 1/2. 
Asymptotic value: fio{oo) — iipo = \/2. 



FIG. 3: (Color online) /io(i) ~ — InP(t) vs t for exponentially 
evanescent traps of a lower initial density po — 0.01. Solid 
line: Eq. 1181 . Symbols: simulation results (there are error 
bars on the symbols but they are too small to see clearly). 
Other parameter values are 7 — 1/2, r — 10*, and K-y = 
l/(2y'7F). Asymptotic value: ^o{(x) = ii/2Po ~ 1.06225. 



we illustrate our earlier caveat, that agreement cannot 
be expected if the initial density of traps is too high 
and the extinction rate of the traps is large, and that 
the agreement improves with lower initial density. The 
disagreement is clear and can be traced exactly to the 
early time trapping events that cumulatively affect the 
survival probability. Fig. |3] shows typical results for the 
lower initial density of traps and a more slowly decay- 
ing trap density, where the agreement between analytic 
results and simulations is clearly very good. 

Finally, it is straightforward to extend the results of 
this section to trap densities that decay as a stretched 
exponential, p{t) ~ po exp[— (i/r)"]. The integral Ijl^ll is 
still straightforward and gives 

an,l2) \' r(,/2a) j' ^^'^ 

which reduces to Eq. H15|) when a = \. The asymptotic 
finite survival probability then is 



P{t) exp 



r(7/2a) 

«r(7/2) 



1 - 



,i-t/rr 



r(7/2a)(t/T)"(i- 



-7/2 

(21) 

An interesting interplay of 7 and a should be noted: there 
are values of a and 7 for which the survival probability 
of the target when the trap density decays as a stretched 
exponential (a < 1) is actually greater than with an ex- 
ponential decay {a — 1). This seemingly counterintuitive 
behavior is connected with the reversal of time inequal- 
ities, i.e., with the fact that {t/r)" is greater (smaller) 
than (t/r) when t is smaller (greater) than r. 



B. Power law decaying trap density 

Suppose now that the trap density decays as a power 
law as might happen, for instance, if there is a process 
of trap-trap annihilation. The trap density at long times 
then decreases as p{t) ~ and it is to be expected 
that the target survival probability (and, in particular, 
whether it is asymptotically vanishing or finite) depends 
sensitively on the relation between the exponents (3 and 
7. We expect that for sufficiently large (3 the target will 
again have a finite probability of surviving forever. 

To find a closed expression for the survival probability 
we need to specify p{t) for all times, not just asymptoti- 
cally, and we choose 



Pit) 



Po 



(l + t/r)'3- 



(22) 



With this form, the integral H15|) can be carried out ex- 
actly, to give 

for all 13, where B is the incomplete Beta function [s^Ib^ 



Ba{z,w) 



with 5R(z) > 0. (24) 



Equation tells us that the typical length explored 
by a (living) trap grows with time as {x'^{t))^/'^ ~ t^/'^. 
On the other hand, the mean distance between traps 
grows as ^ . It thus stands to reason that the 
asymptotic survival probability depends sensitively on 
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the relative magnitudes of (3 and 7/2. To present more 
explicit results in this long-time regime we distinguish 
three cases. 

Case 1: /3 > 7/2. In this case Eq. (|23|) can be written 

as 



^;^i?(7/2,/3-7/2)/_^(7/2,/3-7/2). (25) 



Here B{z,w) is the Beta function (where the requirement 
3fJ(z) > and K(it;) > places us in the "Case 1" regime), 
and Ix{z,w) is the regularized incomplete Beta function 
as defined in Sec. 6.6.2 (pg. 263) of [33. Using the prop- 
erty 6.6.3 in I33 we can set /^(g, b) = 1 — Ii-x{b, a), and 
applying the relation 26.5.5 in [33 we can then write the 
asymptotic result 



/ * (7/2,/3-7/2) = l- 



(t/r) 



7/2-/3 



(/3-7/2)i?(/3-7/2,7/2) 



(26) 

Consequently, recognizing the relation between the Beta 
function and the Gamma function, as t ^ 00 we arrive 
at the asymptotic result 



Mo(^) i'tPO 



r(/3 - 7/2) 
r(/3) 



[t/r) 



7/2-/3 



(/3-7/2)r(7/2) 



(27) 

The survival probability thus approaches (via a power 
law decay of the exponent) the finite asymptotic value 



P{t 00) = exp 



Figure 01 illustrates this result along with numerical sim- 
ulations for comparison. 

Case 2: /3 < 7/2. In this case the integrand in 
Eq. (|13l) goes to zero more slowly than 1/t for t — > 00, 
so that a simple asymptotic analysis of the integral H13|) 
readily establishes that //o(i) goes to infinity with increas- 
ing time as 



Mo(0 



(7/2 - /3)r(7/2) 



7/2-/3 



(29) 



so that the survival probability vanishes at long times as 
a stretched exponential. 



P{t) ^ exp - 



(7/2 - /3)r(7/2) 



{t/r) 



7/2-/3 



(30) 



Analytic and simulation results for this case are shown 
in Fig. |S1 

Case 3: /? = 7/2. This is the marginal case, and 
the incomplete Beta function H24fl can be rewritten as a 
hypergeometric function. 



Mo(i) 



^Po(Vt)^ 
/3r(7/2) 



>i^i(/3,/3,l+/3,-Vr) 



^7P0 

r(7/2) 



Ht/r) 



as t 



(31) 




FIG. 4: (Color online) p,a(t) — — InP(t) vs t for power law 
evanescent traps with /3 > 7/2 ("Case 1"). Parameter values 
are 7 = 0.75, /? = 0.8, po = 0.01, r = 10^ and K-, = 1/(20F). 
Asymptotic value: po{oo) — (.3/4P0 = 2.36549. Solid line: 
Eq. I|28|l . Symbols: simulation results. 
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FIG. 5: (Color online) po{t) — — InP(t) vs t for power law 
evanescent traps with l3 < 7/2 ("Case 2"). Parameter values 
are 7 = 0.8, /3 = 0.2, po = 0.01, r = 10^ and = 1/(20F). 
Solid line: Eq. 1231 . Symbols: simulation results along with 
error bars. 



The survival probability thus decays as an inverse power, 



P{t ^ 00) ^ (t/^)-^Po/r(7/2)_ 



Results for the marginal case are shown in Fig. [S] 



(32) 
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FIG. 6: (Color online) f-ioit) ~ — InP(t) vs t for power law 
evanescent traps with /3 = 7/2 ("Case 3"). Parameter values 
are 7 = 0.8, /? = 0.4, po = 0.01, r = 10^ and = 1/(20F). 
Solid line: Eq. 11231 . Symbols; simulation results. 



IV. CONCLUSIONS 

We have calculated the survival probability of a sta- 
tionary target in a one-dimensional system in which diffu- 
sive or subdiffusive traps that eliminate the target upon 
encounter themselves disappear according to a survival 
probability. The root mean square displacement of the 
traps grows with time as f^^-^, that is, diffusively when 
7 = 1 and subdiffusively when 7 < 1. The survival 
probability of the target depends sensitively on the in- 
terplay of two temporal events, namely, the motion of 
the traps as characterized by the exponent 7 and their 
disappearance. When the motion of the traps is diffu- 
sive or subdiffusive and the traps do not decay in time, 
the survival probability goes to zero as a stretched expo- 
nential, Eq. H14|l . When the traps undergo exponential 
decay or stretched exponential decay, the target has an 
asympototic safety margin, that is, a finite probability of 
surviving forever, cf. Eqs. H19|) and (|21|l . When the traps 
are diffusive or subdiffusive and disappear according to a 
power law survival probability ~ t~^, the survival of the 
target depends sensitively on the relation between 7 and 
(3. If the traps move sufficiently rapidly relative to their 



disappearance, that is, if 7/2 > /?, the target is trapped 
with certainty at long times, its survival probability go- 
ing to zero again as a stretched exponential, cf. Eq. J^OJ). 
If the traps move slowly, 7/2 < /3, then the target has a 
chance of eternal survival, cf. Eq. (|28|l . At the critical 
relation 7/2 = /3 the survival probability goes to zero as 
an inverse power of time, cf. Eq. H32|) . If in fact the trap 
density increases with time, the survival probability of 
the target necessarily vanishes asymptotically. 

In this paper we have calculated the survival proba- 
bility of a target particle in the presence of evanescent 
subdiffusive traps of given time-dependent density. We 
could equally consider the inverse problem, namely, that 
of finding the time dependence of the density of traps 
to obtain a particular survival probability function. For 
this purpose we need only "invert" Eq. (|12|l . 



Pit) 



_ r(7/2) i_^/,p(t) 



Pity 



(33) 



An exponentially decaying survival probability of the 
form P{t) — e~*/^ requires a density that decays as 
p{t) ^ 'Z^. This is included in and consistent with 
Case 2 in Sec. ImBl with 7/2 - /3 = 1. Similarly, for 
an inverse power decay of the form P{t) ^ (t/r)^^ we 
require that p{t) ^ t^^^"^ consistent with Case 3 in the 
same section. 

This work has focused on the survival probability of 
a stationary target. The survival probability of a mov- 
ing target, diffusive or subdiffusive, surrounded by non- 
evanescent diffusive or subdiffusive traps has been consid- 
ered recently in a number of papers pTlll5j| . Extension of 
our work with evanescent traps to the case of a diffusive 
or subdiffusive target is in progress [38l |. 
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